-------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\58592ggr\Arbeit\Projects\Active\VotingIncentives\Data&Ana
> lysis\StataFiles\incetives_in_voting_surveys.smcl
  log type:  smcl
 opened on:   4 Jul 2023, 09:25:29

. 
. gen voted_econ_score = 0 //this uses the party scores of kieskompas

. replace voted_econ_score = voted_econ_score + 1.1 if VVD==1
(280 real changes made)

. replace voted_econ_score = voted_econ_score -0.3 if PVV==1
(279 real changes made)

. replace voted_econ_score = voted_econ_score +0.1 if CDA==1
(125 real changes made)

. replace voted_econ_score = voted_econ_score -0.1 if D66==1
(130 real changes made)

. replace voted_econ_score = voted_econ_score -1.4 if GROEN==1
(137 real changes made)

. replace voted_econ_score = voted_econ_score -2 if SP==1
(165 real changes made)

. replace voted_econ_score = voted_econ_score -1.3 if PvdA==1
(150 real changes made)

. replace voted_econ_score = voted_econ_score - 0.4 if CU==1
(71 real changes made)

. replace voted_econ_score = voted_econ_score -2 if Dieren==1
(125 real changes made)

. replace voted_econ_score = voted_econ_score - 0.4 if PLUS==1
(87 real changes made)

. replace voted_econ_score = voted_econ_score + 0.7 if SGP==1
(17 real changes made)

. replace voted_econ_score = voted_econ_score -1.4 if DENK==1
(16 real changes made)

. replace voted_econ_score = voted_econ_score + 1.4 if Forum==1
(130 real changes made)

. replace voted_econ_score = voted_econ_score - 2 if BIJ1==1
(2 real changes made)

. replace voted_econ_score = voted_econ_score + 0.6 if BBB==1
(3 real changes made)

. 
. gen votes_for_unlisted = Piraten+Oranje+JONG+NIDA+NLBeter+Toekomst+Platform

. gen checked=total_votes - votes_for_unlisted

. 
. replace voted_econ_score = voted_econ_score/(total_votes - votes_for_unlisted
> )
(411 real changes made, 126 to missing)

. 
. gen voted_social_score = 0 //this uses the party scores of kieskompas

. replace voted_social_score = voted_social_score -0.4 if VVD==1
(280 real changes made)

. replace voted_social_score = voted_social_score -1.5 if PVV==1
(279 real changes made)

. replace voted_social_score = voted_social_score -0.5 if CDA==1
(125 real changes made)

. replace voted_social_score = voted_social_score +1.4 if D66==1
(130 real changes made)

. replace voted_social_score = voted_social_score +1.35 if GROEN==1
(137 real changes made)

. replace voted_social_score = voted_social_score +0.8 if SP==1
(165 real changes made)

. replace voted_social_score = voted_social_score +1 if PvdA==1
(150 real changes made)

. replace voted_social_score = voted_social_score +0.3 if CU==1
(71 real changes made)

. replace voted_social_score = voted_social_score +1.3 if Dieren==1
(125 real changes made)

. replace voted_social_score = voted_social_score -0.25 if PLUS==1
(87 real changes made)

. replace voted_social_score = voted_social_score -0.95 if SGP==1
(17 real changes made)

. replace voted_social_score = voted_social_score +0.8 if DENK==1
(16 real changes made)

. replace voted_social_score = voted_social_score -1.5 if Forum==1
(130 real changes made)

. replace voted_social_score = voted_social_score +1.7 if BIJ1==1
(2 real changes made)

. replace voted_social_score = voted_social_score -0.7 if BBB==1
(3 real changes made)

. 
. replace voted_social_score = voted_social_score/(total_votes- votes_for_unlis
> ted)
(411 real changes made, 126 to missing)

. 
. gen ed_center = sqrt(voted_social_score* voted_social_score + voted_econ_scor
> e* voted_econ_score)
(126 missing values generated)

. 
. *scale elicited political views to fit the kieskompas scores
. gen scale_econ=4* econ_views/10 - 2
(124 missing values generated)

. gen scale_social=4* social_views/10 - 2
(124 missing values generated)

. replace scale_social=scale_social*(-1)
(871 real changes made)

. 
. gen diff_econ = abs(scale_econ - voted_econ_score)
(142 missing values generated)

. gen diff_social = abs(scale_social - voted_social_score)
(142 missing values generated)

. 
. gen total_diff = sqrt(diff_econ*diff_econ + diff_social*diff_social)
(142 missing values generated)

. 
. 
. ***Effective number of parties
. set seed 5102020

. 
. capture program drop effective_numberP

. program define effective_numberP, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24
  2. summarize VVD if treatment=="P" & completed==1
  3. scalar `v1' = (r(mean))^2
  4. summarize PVV if treatment=="P" & completed==1
  5. scalar `v2' = (r(mean))^2
  6. summarize CDA if treatment=="P" & completed==1
  7. scalar `v3' = (r(mean))^2
  8. summarize D66 if treatment=="P" & completed==1
  9. scalar `v4' = (r(mean))^2
 10. summarize GROEN if treatment=="P" & completed==1
 11. scalar `v5' = (r(mean))^2
 12. summarize SP if treatment=="P" & completed==1
 13. scalar `v6' = (r(mean))^2
 14. summarize PvdA if treatment=="P" & completed==1
 15. scalar `v7' = (r(mean))^2
 16. summarize CU if treatment=="P" & completed==1
 17. scalar `v8' = (r(mean))^2
 18. summarize Dieren if treatment=="P" & completed==1
 19. scalar `v9' = (r(mean))^2
 20. summarize PLUS if treatment=="P" & completed==1
 21. scalar `v10' = (r(mean))^2
 22. summarize SGP if treatment=="P" & completed==1
 23. scalar `v11' = (r(mean))^2
 24. summarize DENK if treatment=="P" & completed==1
 25. scalar `v12' = (r(mean))^2
 26. summarize Forum if treatment=="P" & completed==1
 27. scalar `v13' = (r(mean))^2
 28. summarize Piraten if treatment=="P" & completed==1
 29. scalar `v14' = (r(mean))^2
 30. summarize BIJ1 if treatment=="P" & completed==1
 31. scalar `v15' = (r(mean))^2
 32. summarize LP if treatment=="P" & completed==1
 33. scalar `v16' = (r(mean))^2
 34. summarize BBB if treatment=="P" & completed==1
 35. scalar `v17' = (r(mean))^2
 36. summarize Oranje if treatment=="P" & completed==1
 37. scalar `v18' = (r(mean))^2
 38. summarize JONG if treatment=="P" & completed==1
 39. scalar `v19' = (r(mean))^2
 40. summarize NIDA if treatment=="P" & completed==1
 41. scalar `v20' = (r(mean))^2
 42. summarize NLBeter if treatment=="P" & completed==1
 43. scalar `v21' = (r(mean))^2
 44. summarize Toekomst if treatment=="P" & completed==1
 45. scalar `v22' = (r(mean))^2
 46. summarize Platform if treatment=="P" & completed==1
 47. scalar `v23' = (r(mean))^2
 48. summarize Volt if treatment=="P" & completed==1
 49. scalar `v24' = (r(mean))^2
 50. return scalar effective_n = 1/(`v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v
> 9'+`v10'+`v11'+`v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v
> 22'+`v23'+`v24')
 51. end

. 
. bootstrap effective_n=r(effective_n) if treatment=="P" & completed==1, reps(1
> 000) nowarn nodots: effective_numberP

Bootstrap results                                        Number of obs =   320
                                                         Replications  = 1,000

      Command: effective_numberP
  effective_n: r(effective_n)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
 effective_n |   9.079624   .4782422    18.99   0.000     8.142287    10.01696
------------------------------------------------------------------------------

. 
. bootstrap effective_n=r(effective_n) if treatment=="P" & completed==1, reps(1
> 000) saving(P_en) nowarn nodots: effective_numberP

Bootstrap results                                        Number of obs =   320
                                                         Replications  = 1,000

      Command: effective_numberP
  effective_n: r(effective_n)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
 effective_n |   9.079624   .4740094    19.15   0.000     8.150583    10.00867
------------------------------------------------------------------------------

. 
. set seed 05102020

. 
. capture program drop effective_numberPI

. program define effective_numberPI, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24
  2. summarize VVD if treatment=="PI" & completed==1
  3. scalar `v1' = (r(mean))^2
  4. summarize PVV if treatment=="PI" & completed==1
  5. scalar `v2' = (r(mean))^2
  6. summarize CDA if treatment=="PI" & completed==1
  7. scalar `v3' = (r(mean))^2
  8. summarize D66 if treatment=="PI" & completed==1
  9. scalar `v4' = (r(mean))^2
 10. summarize GROEN if treatment=="PI" & completed==1
 11. scalar `v5' = (r(mean))^2
 12. summarize SP if treatment=="PI" & completed==1
 13. scalar `v6' = (r(mean))^2
 14. summarize PvdA if treatment=="PI" & completed==1
 15. scalar `v7' = (r(mean))^2
 16. summarize CU if treatment=="PI" & completed==1
 17. scalar `v8' = (r(mean))^2
 18. summarize Dieren if treatment=="PI" & completed==1
 19. scalar `v9' = (r(mean))^2
 20. summarize PLUS if treatment=="PI" & completed==1
 21. scalar `v10' = (r(mean))^2
 22. summarize SGP if treatment=="PI" & completed==1
 23. scalar `v11' = (r(mean))^2
 24. summarize DENK if treatment=="PI" & completed==1
 25. scalar `v12' = (r(mean))^2
 26. summarize Forum if treatment=="PI" & completed==1
 27. scalar `v13' = (r(mean))^2
 28. summarize Piraten if treatment=="PI" & completed==1
 29. scalar `v14' = (r(mean))^2
 30. summarize BIJ1 if treatment=="PI" & completed==1
 31. scalar `v15' = (r(mean))^2
 32. summarize LP if treatment=="PI" & completed==1
 33. scalar `v16' = (r(mean))^2
 34. summarize BBB if treatment=="PI" & completed==1
 35. scalar `v17' = (r(mean))^2
 36. summarize Oranje if treatment=="PI" & completed==1
 37. scalar `v18' = (r(mean))^2
 38. summarize JONG if treatment=="PI" & completed==1
 39. scalar `v19' = (r(mean))^2
 40. summarize NIDA if treatment=="PI" & completed==1
 41. scalar `v20' = (r(mean))^2
 42. summarize NLBeter if treatment=="PI" & completed==1
 43. scalar `v21' = (r(mean))^2
 44. summarize Toekomst if treatment=="PI" & completed==1
 45. scalar `v22' = (r(mean))^2
 46. summarize Platform if treatment=="PI" & completed==1
 47. scalar `v23' = (r(mean))^2
 48. summarize Volt if treatment=="PI" & completed==1
 49. scalar `v24' = (r(mean))^2
 50. return scalar effective_n = 1/(`v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v
> 9'+`v10'+`v11'+`v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v
> 22'+`v23'+`v24')
 51. end

. 
. bootstrap effective_n=r(effective_n) if treatment=="PI" & completed==1, reps(
> 1000) nowarn nodots: effective_numberPI

Bootstrap results                                        Number of obs =   320
                                                         Replications  = 1,000

      Command: effective_numberPI
  effective_n: r(effective_n)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
 effective_n |   10.02742   .5028079    19.94   0.000     9.041933     11.0129
------------------------------------------------------------------------------

. 
. set seed 5102020

. 
. capture program drop effective_numberA

. program define effective_numberA, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24 all
  2. summarize VVD if treatment=="A" & completed==1
  3. scalar `v1' = r(sum)
  4. summarize PVV if treatment=="A" & completed==1
  5. scalar `v2' = r(sum)
  6. summarize CDA if treatment=="A" & completed==1
  7. scalar `v3' = r(sum)
  8. summarize D66 if treatment=="A" & completed==1
  9. scalar `v4' = r(sum)
 10. summarize GROEN if treatment=="A" & completed==1
 11. scalar `v5' = r(sum)
 12. summarize SP if treatment=="A" & completed==1
 13. scalar `v6' = r(sum)
 14. summarize PvdA if treatment=="A" & completed==1
 15. scalar `v7' = r(sum)
 16. summarize CU if treatment=="A" & completed==1
 17. scalar `v8' = r(sum)
 18. summarize Dieren if treatment=="A" & completed==1
 19. scalar `v9' =r(sum)
 20. summarize PLUS if treatment=="A" & completed==1
 21. scalar `v10' = r(sum)
 22. summarize SGP if treatment=="A" & completed==1
 23. scalar `v11' = r(sum)
 24. summarize DENK if treatment=="A" & completed==1
 25. scalar `v12' = r(sum)
 26. summarize Forum if treatment=="A" & completed==1
 27. scalar `v13' = r(sum)
 28. summarize Piraten if treatment=="A" & completed==1
 29. scalar `v14' = r(sum)
 30. summarize BIJ1 if treatment=="A" & completed==1
 31. scalar `v15' = r(sum)
 32. summarize LP if treatment=="A" & completed==1
 33. scalar `v16' = r(sum)
 34. summarize BBB if treatment=="A" & completed==1
 35. scalar `v17' = r(sum)
 36. summarize Oranje if treatment=="A" & completed==1
 37. scalar `v18' = r(sum)
 38. summarize JONG if treatment=="A" & completed==1
 39. scalar `v19' = r(sum)
 40. summarize NIDA if treatment=="A" & completed==1
 41. scalar `v20' = r(sum)
 42. summarize NLBeter if treatment=="A" & completed==1
 43. scalar `v21' = r(sum)
 44. summarize Toekomst if treatment=="A" & completed==1
 45. scalar `v22' = r(sum)
 46. summarize Platform if treatment=="A" & completed==1
 47. scalar `v23' = r(sum)
 48. summarize Volt if treatment=="A" & completed==1
 49. scalar `v24' = r(sum)
 50. scalar `all' = `v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v9'+`v10'+`v11'+`
> v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v22'+`v23'+`v24'
 51. return scalar effective_n = 1/((`v1'/`all')^2+ (`v2'/`all')^2+(`v3'/`all')
> ^2+(`v4'/`all')^2+(`v5'/`all')^2+(`v6'/`all')^2+(`v7'/`all')^2+(`v8'/`all')^2
> +(`v9'/`all')^2+(`v10'/`all')^2+(`v11'/`all')^2+(`v12'/`all')^2+(`v13'/`all')
> ^2+(`v14'/`all')^2+(`v15'/`all')^2+(`v16'/`all')^2+(`v17'/`all')^2+(`v18'/`al
> l')^2+(`v19'/`all')^2+(`v20'/`all')^2+(`v21'/`all')^2+(`v22'/`all')^2+(`v23'/
> `all')^2+(`v24'/`all')^2)
 52. end

. 
. bootstrap effective_n=r(effective_n) if treatment=="A" & completed==1, reps(1
> 000) nowarn nodots: effective_numberA

Bootstrap results                                        Number of obs =   314
                                                         Replications  = 1,000

      Command: effective_numberA
  effective_n: r(effective_n)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
 effective_n |   10.93514   .4422687    24.73   0.000     10.06831    11.80197
------------------------------------------------------------------------------

. 
. set seed 5102020

. 
. capture program drop effective_numberAI

. program define effective_numberAI, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24 all
  2. summarize VVD if treatment=="AI" & completed==1
  3. scalar `v1' = r(sum)
  4. summarize PVV if treatment=="AI" & completed==1
  5. scalar `v2' = r(sum)
  6. summarize CDA if treatment=="AI" & completed==1
  7. scalar `v3' = r(sum)
  8. summarize D66 if treatment=="AI" & completed==1
  9. scalar `v4' = r(sum)
 10. summarize GROEN if treatment=="AI" & completed==1
 11. scalar `v5' = r(sum)
 12. summarize SP if treatment=="AI" & completed==1
 13. scalar `v6' = r(sum)
 14. summarize PvdA if treatment=="AI" & completed==1
 15. scalar `v7' = r(sum)
 16. summarize CU if treatment=="AI" & completed==1
 17. scalar `v8' = r(sum)
 18. summarize Dieren if treatment=="AI" & completed==1
 19. scalar `v9' =r(sum)
 20. summarize PLUS if treatment=="AI" & completed==1
 21. scalar `v10' = r(sum)
 22. summarize SGP if treatment=="AI" & completed==1
 23. scalar `v11' = r(sum)
 24. summarize DENK if treatment=="AI" & completed==1
 25. scalar `v12' = r(sum)
 26. summarize Forum if treatment=="AI" & completed==1
 27. scalar `v13' = r(sum)
 28. summarize Piraten if treatment=="AI" & completed==1
 29. scalar `v14' = r(sum)
 30. summarize BIJ1 if treatment=="AI" & completed==1
 31. scalar `v15' = r(sum)
 32. summarize LP if treatment=="AI" & completed==1
 33. scalar `v16' = r(sum)
 34. summarize BBB if treatment=="AI" & completed==1
 35. scalar `v17' = r(sum)
 36. summarize Oranje if treatment=="AI" & completed==1
 37. scalar `v18' = r(sum)
 38. summarize JONG if treatment=="AI" & completed==1
 39. scalar `v19' = r(sum)
 40. summarize NIDA if treatment=="AI" & completed==1
 41. scalar `v20' = r(sum)
 42. summarize NLBeter if treatment=="AI" & completed==1
 43. scalar `v21' = r(sum)
 44. summarize Toekomst if treatment=="AI" & completed==1
 45. scalar `v22' = r(sum)
 46. summarize Platform if treatment=="AI" & completed==1
 47. scalar `v23' = r(sum)
 48. summarize Volt if treatment=="AI" & completed==1
 49. scalar `v24' = r(sum)
 50. scalar `all' = `v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v9'+`v10'+`v11'+`
> v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v22'+`v23'+`v24'
 51. return scalar effective_n = 1/((`v1'/`all')^2+ (`v2'/`all')^2+(`v3'/`all')
> ^2+(`v4'/`all')^2+(`v5'/`all')^2+(`v6'/`all')^2+(`v7'/`all')^2+(`v8'/`all')^2
> +(`v9'/`all')^2+(`v10'/`all')^2+(`v11'/`all')^2+(`v12'/`all')^2+(`v13'/`all')
> ^2+(`v14'/`all')^2+(`v15'/`all')^2+(`v16'/`all')^2+(`v17'/`all')^2+(`v18'/`al
> l')^2+(`v19'/`all')^2+(`v20'/`all')^2+(`v21'/`all')^2+(`v22'/`all')^2+(`v23'/
> `all')^2+(`v24'/`all')^2)
 52. end

. 
. bootstrap effective_n=r(effective_n) if treatment=="AI" & completed==1, reps(
> 1000) nowarn nodots: effective_numberAI

Bootstrap results                                        Number of obs =   306
                                                         Replications  = 1,000

      Command: effective_numberAI
  effective_n: r(effective_n)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
 effective_n |   11.08681   .4275366    25.93   0.000     10.24885    11.92476
------------------------------------------------------------------------------

. 
. set seed 5102020

. 
. capture program drop differencePPI

. program define differencePPI, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24 vi1 vi2 vi3 vi4 vi5 vi6 vi7 vi8 vi9 vi10 vi11 vi12 vi13
>  vi14 vi15 vi16 vi17 vi18 vi19 vi20 vi21 vi22 vi23 vi24 en eni
  2. summarize VVD if treatment=="P" & completed==1
  3. scalar `v1' = (r(mean))^2
  4. summarize PVV if treatment=="P" & completed==1
  5. scalar `v2' = (r(mean))^2
  6. summarize CDA if treatment=="P" & completed==1
  7. scalar `v3' = (r(mean))^2
  8. summarize D66 if treatment=="P" & completed==1
  9. scalar `v4' = (r(mean))^2
 10. summarize GROEN if treatment=="P" & completed==1
 11. scalar `v5' = (r(mean))^2
 12. summarize SP if treatment=="P" & completed==1
 13. scalar `v6' = (r(mean))^2
 14. summarize PvdA if treatment=="P" & completed==1
 15. scalar `v7' = (r(mean))^2
 16. summarize CU if treatment=="P" & completed==1
 17. scalar `v8' = (r(mean))^2
 18. summarize Dieren if treatment=="P" & completed==1
 19. scalar `v9' = (r(mean))^2
 20. summarize PLUS if treatment=="P" & completed==1
 21. scalar `v10' = (r(mean))^2
 22. summarize SGP if treatment=="P" & completed==1
 23. scalar `v11' = (r(mean))^2
 24. summarize DENK if treatment=="P" & completed==1
 25. scalar `v12' = (r(mean))^2
 26. summarize Forum if treatment=="P" & completed==1
 27. scalar `v13' = (r(mean))^2
 28. summarize Piraten if treatment=="P" & completed==1
 29. scalar `v14' = (r(mean))^2
 30. summarize BIJ1 if treatment=="P" & completed==1
 31. scalar `v15' = (r(mean))^2
 32. summarize LP if treatment=="P" & completed==1
 33. scalar `v16' = (r(mean))^2
 34. summarize BBB if treatment=="P" & completed==1
 35. scalar `v17' = (r(mean))^2
 36. summarize Oranje if treatment=="P" & completed==1
 37. scalar `v18' = (r(mean))^2
 38. summarize JONG if treatment=="P" & completed==1
 39. scalar `v19' = (r(mean))^2
 40. summarize NIDA if treatment=="P" & completed==1
 41. scalar `v20' = (r(mean))^2
 42. summarize NLBeter if treatment=="P" & completed==1
 43. scalar `v21' = (r(mean))^2
 44. summarize Toekomst if treatment=="P" & completed==1
 45. scalar `v22' = (r(mean))^2
 46. summarize Platform if treatment=="P" & completed==1
 47. scalar `v23' = (r(mean))^2
 48. summarize Volt if treatment=="P" & completed==1
 49. scalar `v24' = (r(mean))^2
 50. scalar `en' = 1/(`v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v9'+`v10'+`v11'
> +`v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v22'+`v23'+`v24
> ')
 51. summarize VVD if treatment=="PI" & completed==1
 52. scalar `vi1' = (r(mean))^2
 53. summarize PVV if treatment=="PI" & completed==1
 54. scalar `vi2' = (r(mean))^2
 55. summarize CDA if treatment=="PI" & completed==1
 56. scalar `vi3' = (r(mean))^2
 57. summarize D66 if treatment=="PI" & completed==1
 58. scalar `vi4' = (r(mean))^2
 59. summarize GROEN if treatment=="PI" & completed==1
 60. scalar `vi5' = (r(mean))^2
 61. summarize SP if treatment=="PI" & completed==1
 62. scalar `vi6' = (r(mean))^2
 63. summarize PvdA if treatment=="PI" & completed==1
 64. scalar `vi7' = (r(mean))^2
 65. summarize CU if treatment=="PI" & completed==1
 66. scalar `vi8' = (r(mean))^2
 67. summarize Dieren if treatment=="PI" & completed==1
 68. scalar `vi9' = (r(mean))^2
 69. summarize PLUS if treatment=="PI" & completed==1
 70. scalar `vi10' = (r(mean))^2
 71. summarize SGP if treatment=="PI" & completed==1
 72. scalar `vi11' = (r(mean))^2
 73. summarize DENK if treatment=="PI" & completed==1
 74. scalar `vi12' = (r(mean))^2
 75. summarize Forum if treatment=="PI" & completed==1
 76. scalar `vi13' = (r(mean))^2
 77. summarize Piraten if treatment=="PI" & completed==1
 78. scalar `vi14' = (r(mean))^2
 79. summarize BIJ1 if treatment=="PI" & completed==1
 80. scalar `vi15' = (r(mean))^2
 81. summarize LP if treatment=="PI" & completed==1
 82. scalar `vi16' = (r(mean))^2
 83. summarize BBB if treatment=="PI" & completed==1
 84. scalar `vi17' = (r(mean))^2
 85. summarize Oranje if treatment=="PI" & completed==1
 86. scalar `vi18' = (r(mean))^2
 87. summarize JONG if treatment=="PI" & completed==1
 88. scalar `vi19' = (r(mean))^2
 89. summarize NIDA if treatment=="PI" & completed==1
 90. scalar `vi20' = (r(mean))^2
 91. summarize NLBeter if treatment=="PI" & completed==1
 92. scalar `vi21' = (r(mean))^2
 93. summarize Toekomst if treatment=="PI" & completed==1
 94. scalar `vi22' = (r(mean))^2
 95. summarize Platform if treatment=="PI" & completed==1
 96. scalar `vi23' = (r(mean))^2
 97. summarize Volt if treatment=="PI" & completed==1
 98. scalar `vi24' = (r(mean))^2
 99. scalar `eni' = 1/(`vi1'+ `vi2'+`vi3'+`vi4'+`vi5'+`vi6'+`vi7'+`vi8'+`vi9'+`
> vi10'+`vi11'+`vi12'+`vi13'+`vi14'+`vi15'+`vi16'+`vi17'+`vi18'+`vi19'+`vi20'+`
> vi21'+`vi22'+`vi23'+`vi24')
100. return scalar difference= `en' - `eni'
101. end

. 
. bootstrap difference=r(difference) if treatment=="P" | treatment=="PI", reps(
> 1000) nowarn nodots: differencePPI

Bootstrap results                                        Number of obs =   646
                                                         Replications  = 1,000

      Command: differencePPI
   difference: r(difference)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
  difference |  -.9477947   .7100175    -1.33   0.182    -2.339403     .443814
------------------------------------------------------------------------------

. 
. set seed 5102020

. 
. capture program drop differenceAAI

. program define differenceAAI, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24 vi1 vi2 vi3 vi4 vi5 vi6 vi7 vi8 vi9 vi10 vi11 vi12 vi13
>  vi14 vi15 vi16 vi17 vi18 vi19 vi20 vi21 vi22 vi23 vi24 all en alli eni
  2. summarize VVD if treatment=="A" & completed==1
  3. scalar `v1' = r(sum)
  4. summarize PVV if treatment=="A" & completed==1
  5. scalar `v2' = r(sum)
  6. summarize CDA if treatment=="A" & completed==1
  7. scalar `v3' = r(sum)
  8. summarize D66 if treatment=="A" & completed==1
  9. scalar `v4' = r(sum)
 10. summarize GROEN if treatment=="A" & completed==1
 11. scalar `v5' = r(sum)
 12. summarize SP if treatment=="A" & completed==1
 13. scalar `v6' = r(sum)
 14. summarize PvdA if treatment=="A" & completed==1
 15. scalar `v7' = r(sum)
 16. summarize CU if treatment=="A" & completed==1
 17. scalar `v8' = r(sum)
 18. summarize Dieren if treatment=="A" & completed==1
 19. scalar `v9' =r(sum)
 20. summarize PLUS if treatment=="A" & completed==1
 21. scalar `v10' = r(sum)
 22. summarize SGP if treatment=="A" & completed==1
 23. scalar `v11' = r(sum)
 24. summarize DENK if treatment=="A" & completed==1
 25. scalar `v12' = r(sum)
 26. summarize Forum if treatment=="A" & completed==1
 27. scalar `v13' = r(sum)
 28. summarize Piraten if treatment=="A" & completed==1
 29. scalar `v14' = r(sum)
 30. summarize BIJ1 if treatment=="A" & completed==1
 31. scalar `v15' = r(sum)
 32. summarize LP if treatment=="A" & completed==1
 33. scalar `v16' = r(sum)
 34. summarize BBB if treatment=="A" & completed==1
 35. scalar `v17' = r(sum)
 36. summarize Oranje if treatment=="A" & completed==1
 37. scalar `v18' = r(sum)
 38. summarize JONG if treatment=="A" & completed==1
 39. scalar `v19' = r(sum)
 40. summarize NIDA if treatment=="A" & completed==1
 41. scalar `v20' = r(sum)
 42. summarize NLBeter if treatment=="A" & completed==1
 43. scalar `v21' = r(sum)
 44. summarize Toekomst if treatment=="A" & completed==1
 45. scalar `v22' = r(sum)
 46. summarize Platform if treatment=="A" & completed==1
 47. scalar `v23' = r(sum)
 48. summarize Volt if treatment=="A" & completed==1
 49. scalar `v24' = r(sum)
 50. scalar `all' = `v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v9'+`v10'+`v11'+`
> v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v22'+`v23'+`v24'
 51. scalar `en' = 1/((`v1'/`all')^2+ (`v2'/`all')^2+(`v3'/`all')^2+(`v4'/`all'
> )^2+(`v5'/`all')^2+(`v6'/`all')^2+(`v7'/`all')^2+(`v8'/`all')^2+(`v9'/`all')^
> 2+(`v10'/`all')^2+(`v11'/`all')^2+(`v12'/`all')^2+(`v13'/`all')^2+(`v14'/`all
> ')^2+(`v15'/`all')^2+(`v16'/`all')^2+(`v17'/`all')^2+(`v18'/`all')^2+(`v19'/`
> all')^2+(`v20'/`all')^2+(`v21'/`all')^2+(`v22'/`all')^2+(`v23'/`all')^2+(`v24
> '/`all')^2)
 52. summarize VVD if treatment=="AI" & completed==1
 53. scalar `vi1' = r(sum)
 54. summarize PVV if treatment=="AI" & completed==1
 55. scalar `vi2' = r(sum)
 56. summarize CDA if treatment=="AI" & completed==1
 57. scalar `vi3' = r(sum)
 58. summarize D66 if treatment=="AI" & completed==1
 59. scalar `vi4' = r(sum)
 60. summarize GROEN if treatment=="AI" & completed==1
 61. scalar `vi5' = r(sum)
 62. summarize SP if treatment=="AI" & completed==1
 63. scalar `vi6' = r(sum)
 64. summarize PvdA if treatment=="AI" & completed==1
 65. scalar `vi7' = r(sum)
 66. summarize CU if treatment=="AI" & completed==1
 67. scalar `vi8' = r(sum)
 68. summarize Dieren if treatment=="AI" & completed==1
 69. scalar `vi9' =r(sum)
 70. summarize PLUS if treatment=="AI" & completed==1
 71. scalar `vi10' = r(sum)
 72. summarize SGP if treatment=="AI" & completed==1
 73. scalar `vi11' = r(sum)
 74. summarize DENK if treatment=="AI" & completed==1
 75. scalar `vi12' = r(sum)
 76. summarize Forum if treatment=="AI" & completed==1
 77. scalar `vi13' = r(sum)
 78. summarize Piraten if treatment=="AI" & completed==1
 79. scalar `vi14' = r(sum)
 80. summarize BIJ1 if treatment=="AI" & completed==1
 81. scalar `vi15' = r(sum)
 82. summarize LP if treatment=="AI" & completed==1
 83. scalar `vi16' = r(sum)
 84. summarize BBB if treatment=="AI" & completed==1
 85. scalar `vi17' = r(sum)
 86. summarize Oranje if treatment=="AI" & completed==1
 87. scalar `vi18' = r(sum)
 88. summarize JONG if treatment=="AI" & completed==1
 89. scalar `vi19' = r(sum)
 90. summarize NIDA if treatment=="AI" & completed==1
 91. scalar `vi20' = r(sum)
 92. summarize NLBeter if treatment=="AI" & completed==1
 93. scalar `vi21' = r(sum)
 94. summarize Toekomst if treatment=="AI" & completed==1
 95. scalar `vi22' = r(sum)
 96. summarize Platform if treatment=="AI" & completed==1
 97. scalar `vi23' = r(sum)
 98. summarize Volt if treatment=="AI" & completed==1
 99. scalar `vi24' = r(sum)
100. scalar `alli' = `vi1'+ `vi2'+`vi3'+`vi4'+`vi5'+`vi6'+`vi7'+`vi8'+`vi9'+`vi
> 10'+`vi11'+`vi12'+`vi13'+`vi14'+`vi15'+`vi16'+`vi17'+`vi18'+`vi19'+`vi20'+`vi
> 21'+`vi22'+`vi23'+`vi24'
101. scalar `eni' = 1/((`vi1'/`alli')^2+ (`vi2'/`alli')^2+(`vi3'/`alli')^2+(`vi
> 4'/`alli')^2+(`vi5'/`alli')^2+(`vi6'/`alli')^2+(`vi7'/`alli')^2+(`vi8'/`alli'
> )^2+(`vi9'/`alli')^2+(`vi10'/`alli')^2+(`vi11'/`alli')^2+(`vi12'/`alli')^2+(`
> vi13'/`alli')^2+(`vi14'/`alli')^2+(`vi15'/`alli')^2+(`vi16'/`alli')^2+(`vi17'
> /`alli')^2+(`vi18'/`alli')^2+(`vi19'/`alli')^2+(`vi20'/`alli')^2+(`vi21'/`all
> i')^2+(`vi22'/`alli')^2+(`vi23'/`alli')^2+(`vi24'/`alli')^2)
102. return scalar difference = `en'-`eni'
103. end

. 
. bootstrap difference=r(difference) if treatment=="A" | treatment=="AI", reps(
> 1000) nowarn nodots: differenceAAI

Bootstrap results                                        Number of obs =   633
                                                         Replications  = 1,000

      Command: differenceAAI
   difference: r(difference)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
  difference |  -.1516672   .6192846    -0.24   0.807    -1.365443    1.062108
------------------------------------------------------------------------------

. 
. set seed 5102020

. 
. capture program drop differencePA

. program define differencePA, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24 vi1 vi2 vi3 vi4 vi5 vi6 vi7 vi8 vi9 vi10 vi11 vi12 vi13
>  vi14 vi15 vi16 vi17 vi18 vi19 vi20 vi21 vi22 vi23 vi24 en eni alli
  2. summarize VVD if treatment_num==1 & completed==1
  3. scalar `v1' = (r(mean))^2
  4. summarize PVV if treatment_num==1 & completed==1
  5. scalar `v2' = (r(mean))^2
  6. summarize CDA if treatment_num==1 & completed==1
  7. scalar `v3' = (r(mean))^2
  8. summarize D66 if treatment_num==1 & completed==1
  9. scalar `v4' = (r(mean))^2
 10. summarize GROEN if treatment_num==1 & completed==1
 11. scalar `v5' = (r(mean))^2
 12. summarize SP if treatment_num==1 & completed==1
 13. scalar `v6' = (r(mean))^2
 14. summarize PvdA if treatment_num==1 & completed==1
 15. scalar `v7' = (r(mean))^2
 16. summarize CU if treatment_num==1 & completed==1
 17. scalar `v8' = (r(mean))^2
 18. summarize Dieren if treatment_num==1 & completed==1
 19. scalar `v9' = (r(mean))^2
 20. summarize PLUS if treatment_num==1 & completed==1
 21. scalar `v10' = (r(mean))^2
 22. summarize SGP if treatment_num==1 & completed==1
 23. scalar `v11' = (r(mean))^2
 24. summarize DENK if treatment_num==1 & completed==1
 25. scalar `v12' = (r(mean))^2
 26. summarize Forum if treatment_num==1 & completed==1
 27. scalar `v13' = (r(mean))^2
 28. summarize Piraten if treatment_num==1 & completed==1
 29. scalar `v14' = (r(mean))^2
 30. summarize BIJ1 if treatment_num==1 & completed==1
 31. scalar `v15' = (r(mean))^2
 32. summarize LP if treatment_num==1 & completed==1
 33. scalar `v16' = (r(mean))^2
 34. summarize BBB if treatment_num==1 & completed==1
 35. scalar `v17' = (r(mean))^2
 36. summarize Oranje if treatment_num==1 & completed==1
 37. scalar `v18' = (r(mean))^2
 38. summarize JONG if treatment_num==1 & completed==1
 39. scalar `v19' = (r(mean))^2
 40. summarize NIDA if treatment_num==1 & completed==1
 41. scalar `v20' = (r(mean))^2
 42. summarize NLBeter if treatment_num==1 & completed==1
 43. scalar `v21' = (r(mean))^2
 44. summarize Toekomst if treatment_num==1 & completed==1
 45. scalar `v22' = (r(mean))^2
 46. summarize Platform if treatment_num==1 & completed==1
 47. scalar `v23' = (r(mean))^2
 48. summarize Volt if treatment_num==1 & completed==1
 49. scalar `v24' = (r(mean))^2
 50. scalar `en' = 1/(`v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v9'+`v10'+`v11'
> +`v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v22'+`v23'+`v24
> ')
 51. summarize VVD if treatment_num==3 & completed==1
 52. scalar `vi1' = r(sum)
 53. summarize PVV if treatment_num==3 & completed==1
 54. scalar `vi2' = r(sum)
 55. summarize CDA if treatment_num==3 & completed==1
 56. scalar `vi3' = r(sum)
 57. summarize D66 if treatment_num==3 & completed==1
 58. scalar `vi4' = r(sum)
 59. summarize GROEN if treatment_num==3 & completed==1
 60. scalar `vi5' = r(sum)
 61. summarize SP if treatment_num==3 & completed==1
 62. scalar `vi6' = r(sum)
 63. summarize PvdA if treatment_num==3 & completed==1
 64. scalar `vi7' = r(sum)
 65. summarize CU if treatment_num==3 & completed==1
 66. scalar `vi8' = r(sum)
 67. summarize Dieren if treatment_num==3 & completed==1
 68. scalar `vi9' =r(sum)
 69. summarize PLUS if treatment_num==3 & completed==1
 70. scalar `vi10' = r(sum)
 71. summarize SGP if treatment_num==3 & completed==1
 72. scalar `vi11' = r(sum)
 73. summarize DENK if treatment_num==3 & completed==1
 74. scalar `vi12' = r(sum)
 75. summarize Forum if treatment_num==3 & completed==1
 76. scalar `vi13' = r(sum)
 77. summarize Piraten if treatment_num==3 & completed==1
 78. scalar `vi14' = r(sum)
 79. summarize BIJ1 if treatment_num==3 & completed==1
 80. scalar `vi15' = r(sum)
 81. summarize LP if treatment_num==3 & completed==1
 82. scalar `vi16' = r(sum)
 83. summarize BBB if treatment_num==3 & completed==1
 84. scalar `vi17' = r(sum)
 85. summarize Oranje if treatment_num==3 & completed==1
 86. scalar `vi18' = r(sum)
 87. summarize JONG if treatment_num==3 & completed==1
 88. scalar `vi19' = r(sum)
 89. summarize NIDA if treatment_num==3 & completed==1
 90. scalar `vi20' = r(sum)
 91. summarize NLBeter if treatment_num==3 & completed==1
 92. scalar `vi21' = r(sum)
 93. summarize Toekomst if treatment_num==3 & completed==1
 94. scalar `vi22' = r(sum)
 95. summarize Platform if treatment_num==3 & completed==1
 96. scalar `vi23' = r(sum)
 97. summarize Volt if treatment_num==3 & completed==1
 98. scalar `vi24' = r(sum)
 99. scalar `alli' = `vi1'+ `vi2'+`vi3'+`vi4'+`vi5'+`vi6'+`vi7'+`vi8'+`vi9'+`vi
> 10'+`vi11'+`vi12'+`vi13'+`vi14'+`vi15'+`vi16'+`vi17'+`vi18'+`vi19'+`vi20'+`vi
> 21'+`vi22'+`vi23'+`vi24'
100. scalar `eni' = 1/((`vi1'/`alli')^2+ (`vi2'/`alli')^2+(`vi3'/`alli')^2+(`vi
> 4'/`alli')^2+(`vi5'/`alli')^2+(`vi6'/`alli')^2+(`vi7'/`alli')^2+(`vi8'/`alli'
> )^2+(`vi9'/`alli')^2+(`vi10'/`alli')^2+(`vi11'/`alli')^2+(`vi12'/`alli')^2+(`
> vi13'/`alli')^2+(`vi14'/`alli')^2+(`vi15'/`alli')^2+(`vi16'/`alli')^2+(`vi17'
> /`alli')^2+(`vi18'/`alli')^2+(`vi19'/`alli')^2+(`vi20'/`alli')^2+(`vi21'/`all
> i')^2+(`vi22'/`alli')^2+(`vi23'/`alli')^2+(`vi24'/`alli')^2)
101. return scalar difference= `en' - `eni'
102. end

. 
. **permute treatment_num r(difference) , reps(1000) dots(100): differencePA if
>  treatment_num==1 | treatment_num==3
. 
. bootstrap difference=r(difference) if treatment=="P" | treatment=="A", reps(1
> 000) nowarn nodots: differencePA

Bootstrap results                                        Number of obs =   640
                                                         Replications  = 1,000

      Command: differencePA
   difference: r(difference)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
  difference |  -1.855515    .645747    -2.87   0.004    -3.121156   -.5898746
------------------------------------------------------------------------------

. 
. set seed 5102020

. 
. capture program drop differencePIAI

. program define differencePIAI, rclass
  1. tempname v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17 v18 v1
> 9 v20 v21 v22 v23 v24 vi1 vi2 vi3 vi4 vi5 vi6 vi7 vi8 vi9 vi10 vi11 vi12 vi13
>  vi14 vi15 vi16 vi17 vi18 vi19 vi20 vi21 vi22 vi23 vi24 en eni alli
  2. summarize VVD if treatment=="PI" & completed==1
  3. scalar `v1' = (r(mean))^2
  4. summarize PVV if treatment=="PI" & completed==1
  5. scalar `v2' = (r(mean))^2
  6. summarize CDA if treatment=="PI" & completed==1
  7. scalar `v3' = (r(mean))^2
  8. summarize D66 if treatment=="PI" & completed==1
  9. scalar `v4' = (r(mean))^2
 10. summarize GROEN if treatment=="PI" & completed==1
 11. scalar `v5' = (r(mean))^2
 12. summarize SP if treatment=="PI" & completed==1
 13. scalar `v6' = (r(mean))^2
 14. summarize PvdA if treatment=="PI" & completed==1
 15. scalar `v7' = (r(mean))^2
 16. summarize CU if treatment=="PI" & completed==1
 17. scalar `v8' = (r(mean))^2
 18. summarize Dieren if treatment=="PI" & completed==1
 19. scalar `v9' = (r(mean))^2
 20. summarize PLUS if treatment=="PI" & completed==1
 21. scalar `v10' = (r(mean))^2
 22. summarize SGP if treatment=="PI" & completed==1
 23. scalar `v11' = (r(mean))^2
 24. summarize DENK if treatment=="PI" & completed==1
 25. scalar `v12' = (r(mean))^2
 26. summarize Forum if treatment=="PI" & completed==1
 27. scalar `v13' = (r(mean))^2
 28. summarize Piraten if treatment=="PI" & completed==1
 29. scalar `v14' = (r(mean))^2
 30. summarize BIJ1 if treatment=="PI" & completed==1
 31. scalar `v15' = (r(mean))^2
 32. summarize LP if treatment=="PI" & completed==1
 33. scalar `v16' = (r(mean))^2
 34. summarize BBB if treatment=="PI" & completed==1
 35. scalar `v17' = (r(mean))^2
 36. summarize Oranje if treatment=="PI" & completed==1
 37. scalar `v18' = (r(mean))^2
 38. summarize JONG if treatment=="PI" & completed==1
 39. scalar `v19' = (r(mean))^2
 40. summarize NIDA if treatment=="PI" & completed==1
 41. scalar `v20' = (r(mean))^2
 42. summarize NLBeter if treatment=="PI" & completed==1
 43. scalar `v21' = (r(mean))^2
 44. summarize Toekomst if treatment=="PI" & completed==1
 45. scalar `v22' = (r(mean))^2
 46. summarize Platform if treatment=="PI" & completed==1
 47. scalar `v23' = (r(mean))^2
 48. summarize Volt if treatment=="PI" & completed==1
 49. scalar `v24' = (r(mean))^2
 50. scalar `en' = 1/(`v1'+ `v2'+`v3'+`v4'+`v5'+`v6'+`v7'+`v8'+`v9'+`v10'+`v11'
> +`v12'+`v13'+`v14'+`v15'+`v16'+`v17'+`v18'+`v19'+`v20'+`v21'+`v22'+`v23'+`v24
> ')
 51. summarize VVD if treatment=="AI" & completed==1
 52. scalar `vi1' = r(sum)
 53. summarize PVV if treatment=="AI" & completed==1
 54. scalar `vi2' = r(sum)
 55. summarize CDA if treatment=="AI" & completed==1
 56. scalar `vi3' = r(sum)
 57. summarize D66 if treatment=="AI" & completed==1
 58. scalar `vi4' = r(sum)
 59. summarize GROEN if treatment=="AI" & completed==1
 60. scalar `vi5' = r(sum)
 61. summarize SP if treatment=="AI" & completed==1
 62. scalar `vi6' = r(sum)
 63. summarize PvdA if treatment=="AI" & completed==1
 64. scalar `vi7' = r(sum)
 65. summarize CU if treatment=="AI" & completed==1
 66. scalar `vi8' = r(sum)
 67. summarize Dieren if treatment=="AI" & completed==1
 68. scalar `vi9' =r(sum)
 69. summarize PLUS if treatment=="AI" & completed==1
 70. scalar `vi10' = r(sum)
 71. summarize SGP if treatment=="AI" & completed==1
 72. scalar `vi11' = r(sum)
 73. summarize DENK if treatment=="AI" & completed==1
 74. scalar `vi12' = r(sum)
 75. summarize Forum if treatment=="AI" & completed==1
 76. scalar `vi13' = r(sum)
 77. summarize Piraten if treatment=="AI" & completed==1
 78. scalar `vi14' = r(sum)
 79. summarize BIJ1 if treatment=="AI" & completed==1
 80. scalar `vi15' = r(sum)
 81. summarize LP if treatment=="AI" & completed==1
 82. scalar `vi16' = r(sum)
 83. summarize BBB if treatment=="AI" & completed==1
 84. scalar `vi17' = r(sum)
 85. summarize Oranje if treatment=="AI" & completed==1
 86. scalar `vi18' = r(sum)
 87. summarize JONG if treatment=="AI" & completed==1
 88. scalar `vi19' = r(sum)
 89. summarize NIDA if treatment=="AI" & completed==1
 90. scalar `vi20' = r(sum)
 91. summarize NLBeter if treatment=="AI" & completed==1
 92. scalar `vi21' = r(sum)
 93. summarize Toekomst if treatment=="AI" & completed==1
 94. scalar `vi22' = r(sum)
 95. summarize Platform if treatment=="AI" & completed==1
 96. scalar `vi23' = r(sum)
 97. summarize Volt if treatment=="AI" & completed==1
 98. scalar `vi24' = r(sum)
 99. scalar `alli' = `vi1'+ `vi2'+`vi3'+`vi4'+`vi5'+`vi6'+`vi7'+`vi8'+`vi9'+`vi
> 10'+`vi11'+`vi12'+`vi13'+`vi14'+`vi15'+`vi16'+`vi17'+`vi18'+`vi19'+`vi20'+`vi
> 21'+`vi22'+`vi23'+`vi24'
100. scalar `eni' = 1/((`vi1'/`alli')^2+ (`vi2'/`alli')^2+(`vi3'/`alli')^2+(`vi
> 4'/`alli')^2+(`vi5'/`alli')^2+(`vi6'/`alli')^2+(`vi7'/`alli')^2+(`vi8'/`alli'
> )^2+(`vi9'/`alli')^2+(`vi10'/`alli')^2+(`vi11'/`alli')^2+(`vi12'/`alli')^2+(`
> vi13'/`alli')^2+(`vi14'/`alli')^2+(`vi15'/`alli')^2+(`vi16'/`alli')^2+(`vi17'
> /`alli')^2+(`vi18'/`alli')^2+(`vi19'/`alli')^2+(`vi20'/`alli')^2+(`vi21'/`all
> i')^2+(`vi22'/`alli')^2+(`vi23'/`alli')^2+(`vi24'/`alli')^2)
101. return scalar difference= `en' - `eni'
102. end

. 
. bootstrap difference=r(difference), reps(1000) nowarn nodots: differencePIAI

Bootstrap results                                        Number of obs = 1,386
                                                         Replications  = 1,000

      Command: differencePIAI
   difference: r(difference)

------------------------------------------------------------------------------
             |   Observed   Bootstrap                         Normal-based
             | coefficient  std. err.      z    P>|z|     [95% conf. interval]
-------------+----------------------------------------------------------------
  difference |  -1.059388   .6418647    -1.65   0.099     -2.31742    .1986437
------------------------------------------------------------------------------

. 
. tabi 60 53 20 20 26 34 33 8 24 9 5 0 19 0 0 0 1 0 0 0 0 6 0 2 \ 51 51 20 22 1
> 7 32 26 12 32 16 3 2 26 2 0 1 0 1 1 0 0 4 0 1

           |                                                                   
>               col
       row |         1          2          3          4          5          6  
>         7          8          9         10         11         12         13  
>        14         16 |     Total
-----------+-------------------------------------------------------------------
-----------------------------------------------------------------------------
---------------------+----------
         1 |        60         53         20         20         26         34  
>        33          8         24          9          5          0         19  
>         0          0 |       320 
         2 |        51         51         20         22         17         32  
>        26         12         32         16          3          2         26  
>         2          1 |       320 
-----------+-------------------------------------------------------------------
-----------------------------------------------------------------------------
---------------------+----------
     Total |       111        104         40         42         43         66  
>        59         20         56         25          8          2         45  
>         2          1 |       640 


           |                          col
       row |        17         18         19         22         24 |     Total
-----------+-------------------------------------------------------+----------
         1 |         1          0          0          6          2 |       320 
         2 |         0          1          1          4          1 |       320 
-----------+-------------------------------------------------------+----------
     Total |         1          1          1         10          3 |       640 

         Pearson chi2(19) =  17.8633   Pr = 0.532

. tabi 95 89 39 45 49 45 43 24 38 32 7 10 50 4 1 2 1 3 1 1 2 8 4 4 \ 73 83 42 4
> 0 42 49 45 27 29 29 2 4 34 5 1 2 1 2 2 1 3 10 4 2

           |                                                                   
>               col
       row |         1          2          3          4          5          6  
>         7          8          9         10         11         12         13  
>        14         15 |     Total
-----------+-------------------------------------------------------------------
-----------------------------------------------------------------------------
---------------------+----------
         1 |        95         89         39         45         49         45  
>        43         24         38         32          7         10         50  
>         4          1 |       597 
         2 |        73         83         42         40         42         49  
>        45         27         29         29          2          4         34  
>         5          1 |       532 
-----------+-------------------------------------------------------------------
-----------------------------------------------------------------------------
---------------------+----------
     Total |       168        172         81         85         91         94  
>        88         51         67         61          9         14         84  
>         9          2 |     1,129 


           |                                                col
       row |        16         17         18         19         20         21  
>        22         23         24 |     Total
-----------+-------------------------------------------------------------------
--------------------------------+----------
         1 |         2          1          3          1          1          2  
>         8          4          4 |       597 
         2 |         2          1          2          2          1          3  
>        10          4          2 |       532 
-----------+-------------------------------------------------------------------
--------------------------------+----------
     Total |         4          2          5          3          2          5  
>        18          8          6 |     1,129 

         Pearson chi2(23) =  12.2110   Pr = 0.967

. ranksum total_votes if completed==1 & (treatment=="A" | treatment=="AI"), by(
> treatment)

Two-sample Wilcoxon rank-sum (Mann–Whitney) test

   treatment |      Obs    Rank sum    Expected
-------------+---------------------------------
           A |      314    101163.5       97497
          AI |      306     91346.5       95013
-------------+---------------------------------
    Combined |      620      192510      192510

Unadjusted variance  4972347.00
Adjustment for ties  -842378.75
                     ----------
Adjusted variance    4129968.25

H0: total_~s(treatm~t==A) = total_~s(treatm~t==AI)
         z =  1.804
Prob > |z| = 0.0712

Note: Exact p-value is not computed by default for sample sizes > 200.
      Use option exact to compute it.

. ttest total_votes if completed==1 & (treatment=="A" | treatment=="AI"), by(tr
> eatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       A |     314    1.901274    .0675107    1.196293    1.768442    2.034106
      AI |     306    1.738562    .0611582    1.069832    1.618217    1.858908
---------+--------------------------------------------------------------------
Combined |     620    1.820968    .0456888    1.137642    1.731244    1.910692
---------+--------------------------------------------------------------------
    diff |            .1627118    .0912246                -.016436    .3418595
------------------------------------------------------------------------------
    diff = mean(A) - mean(AI)                                     t =   1.7836
H0: diff = 0                                     Degrees of freedom =      618

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.9625         Pr(|T| > |t|) = 0.0750          Pr(T > t) = 0.0375

. 
. ***Analyse centrality of votes
. tabulate treatment if completed==1 & checked!=0 , sum( ed_center)

  Treatment |
 as string: |
  P, PI, A, |        Summary of ed_center
         AI |        Mean   Std. dev.       Freq.
------------+------------------------------------
          A |   1.2472759   .57584891         311
         AI |   1.2702434   .57110523         305
          P |   1.5248467   .54789907         314
         PI |   1.5173498   .59368787         312
------------+------------------------------------
      Total |   1.3909356   .58658433       1,242

. ci mean ed_center if completed==1 & votes_for_unlisted==0 & treatment=="P"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
   ed_center |        314    1.524847    .0309197         1.46401    1.585684

. ci mean ed_center if completed==1 & votes_for_unlisted==0 & treatment=="PI"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
   ed_center |        312     1.51735     .033611        1.451216    1.583483

. ci mean ed_center if completed==1 & votes_for_unlisted==0 & treatment=="A"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
   ed_center |        296    1.252971    .0334968        1.187048    1.318894

. ci mean ed_center if completed==1 & votes_for_unlisted==0 & treatment=="AI"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
   ed_center |        282    1.296978    .0332561        1.231515     1.36244

. ttest ed_center if completed==1 & votes_for_unlisted==0 & (treatment=="P" | t
> reatment=="PI"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       P |     314    1.524847    .0309197    .5478991     1.46401    1.585684
      PI |     312     1.51735     .033611    .5936879    1.451216    1.583483
---------+--------------------------------------------------------------------
Combined |     626     1.52111    .0228111    .5707343    1.476314    1.565906
---------+--------------------------------------------------------------------
    diff |             .007497     .045658               -.0821651     .097159
------------------------------------------------------------------------------
    diff = mean(P) - mean(PI)                                     t =   0.1642
H0: diff = 0                                     Degrees of freedom =      624

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.5652         Pr(|T| > |t|) = 0.8696          Pr(T > t) = 0.4348

. ttest ed_center if completed==1 & votes_for_unlisted==0 & (treatment=="A" | t
> reatment=="AI"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       A |     296    1.252971    .0334968    .5763007    1.187048    1.318894
      AI |     282    1.296978    .0332561    .5584652    1.231515     1.36244
---------+--------------------------------------------------------------------
Combined |     578    1.274441    .0236092    .5676047    1.228071    1.320812
---------+--------------------------------------------------------------------
    diff |           -.0440064    .0472377               -.1367856    .0487728
------------------------------------------------------------------------------
    diff = mean(A) - mean(AI)                                     t =  -0.9316
H0: diff = 0                                     Degrees of freedom =      576

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1760         Pr(|T| > |t|) = 0.3519          Pr(T > t) = 0.8240

. ttest ed_center if completed==1 & votes_for_unlisted==0 & (treatment=="P" | t
> reatment=="A"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       A |     296    1.252971    .0334968    .5763007    1.187048    1.318894
       P |     314    1.524847    .0309197    .5478991     1.46401    1.585684
---------+--------------------------------------------------------------------
Combined |     610     1.39292    .0233877    .5776333     1.34699    1.438851
---------+--------------------------------------------------------------------
    diff |           -.2718755    .0455178               -.3612667   -.1824843
------------------------------------------------------------------------------
    diff = mean(A) - mean(P)                                      t =  -5.9729
H0: diff = 0                                     Degrees of freedom =      608

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0000         Pr(|T| > |t|) = 0.0000          Pr(T > t) = 1.0000

. ttest ed_center if completed==1 & votes_for_unlisted==0 & (treatment=="PI" | 
> treatment=="AI"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
Variable |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
      AI |     282    1.296978    .0332561    .5584652    1.231515     1.36244
      PI |     312     1.51735     .033611    .5936879    1.451216    1.583483
---------+--------------------------------------------------------------------
Combined |     594    1.412729     .024092    .5871721    1.365413    1.460045
---------+--------------------------------------------------------------------
    diff |           -.2203722    .0474292               -.3135221   -.1272222
------------------------------------------------------------------------------
    diff = mean(AI) - mean(PI)                                    t =  -4.6463
H0: diff = 0                                     Degrees of freedom =      592

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0000         Pr(|T| > |t|) = 0.0000          Pr(T > t) = 1.0000

. 
. tabulate treatment if completed==1 & votes_for_unlisted==0 , sum( total_diff)

  Treatment |
 as string: |
  P, PI, A, |        Summary of total_diff
         AI |        Mean   Std. dev.       Freq.
------------+------------------------------------
          A |   1.3566986    .7630936         296
         AI |   1.4212566   .81067935         282
          P |   1.4796369    .8754585         314
         PI |   1.5183657    .8283065         312
------------+------------------------------------
      Total |   1.4457751   .82270567       1,204

. ci mean total_diff if completed==1 & votes_for_unlisted==0 & treatment=="P"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
  total_diff |        314    1.479637     .049405        1.382429    1.576845

. ci mean total_diff if completed==1 & votes_for_unlisted==0 & treatment=="PI"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
  total_diff |        312    1.518366    .0468936        1.426097    1.610635

. ci mean total_diff if completed==1 & votes_for_unlisted==0 & treatment=="A"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
  total_diff |        296    1.356699    .0443539        1.269408    1.443989

. ci mean total_diff if completed==1 & votes_for_unlisted==0 & treatment=="AI"

    Variable |        Obs        Mean    Std. err.       [95% conf. interval]
-------------+---------------------------------------------------------------
  total_diff |        282    1.421257    .0482753         1.32623    1.516284

. ttest total_diff if completed==1 & votes_for_unlisted==0 & (treatment=="P" | 
> treatment=="PI"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       P |     314    1.479637     .049405    .8754585    1.382429    1.576845
      PI |     312    1.518366    .0468936    .8283065    1.426097    1.610635
---------+--------------------------------------------------------------------
Combined |     626    1.498939    .0340457    .8518226    1.432082    1.565797
---------+--------------------------------------------------------------------
    diff |           -.0387288    .0681286               -.1725179    .0950604
------------------------------------------------------------------------------
    diff = mean(P) - mean(PI)                                     t =  -0.5685
H0: diff = 0                                     Degrees of freedom =      624

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.2850         Pr(|T| > |t|) = 0.5699          Pr(T > t) = 0.7150

. ttest total_diff if completed==1 & votes_for_unlisted==0 & (treatment=="A" | 
> treatment=="AI"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       A |     296    1.356699    .0443539    .7630936    1.269408    1.443989
      AI |     282    1.421257    .0482753    .8106793     1.32623    1.516284
---------+--------------------------------------------------------------------
Combined |     578    1.388196    .0327203    .7866492     1.32393    1.452461
---------+--------------------------------------------------------------------
    diff |            -.064558    .0654613               -.1931301     .064014
------------------------------------------------------------------------------
    diff = mean(A) - mean(AI)                                     t =  -0.9862
H0: diff = 0                                     Degrees of freedom =      576

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.1622         Pr(|T| > |t|) = 0.3244          Pr(T > t) = 0.8378

. ttest total_diff if completed==1 & votes_for_unlisted==0 & (treatment=="P" | 
> treatment=="A"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
       A |     296    1.356699    .0443539    .7630936    1.269408    1.443989
       P |     314    1.479637     .049405    .8754585    1.382429    1.576845
---------+--------------------------------------------------------------------
Combined |     610    1.419982    .0333822    .8244784    1.354423     1.48554
---------+--------------------------------------------------------------------
    diff |           -.1229384    .0666621               -.2538543    .0079776
------------------------------------------------------------------------------
    diff = mean(A) - mean(P)                                      t =  -1.8442
H0: diff = 0                                     Degrees of freedom =      608

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0328         Pr(|T| > |t|) = 0.0656          Pr(T > t) = 0.9672

. ttest total_diff if completed==1 & votes_for_unlisted==0 & (treatment=="PI" |
>  treatment=="AI"), by(treatment)

Two-sample t test with equal variances
------------------------------------------------------------------------------
Variable |     Obs        Mean    Std. err.   Std. dev.   [95% conf. interval]
---------+--------------------------------------------------------------------
      AI |     282    1.421257    .0482753    .8106793     1.32623    1.516284
      PI |     312    1.518366    .0468936    .8283065    1.426097    1.610635
---------+--------------------------------------------------------------------
Combined |     594    1.472263     .033675    .8207314    1.406127      1.5384
---------+--------------------------------------------------------------------
    diff |           -.0971091    .0673749                -.229432    .0352139
------------------------------------------------------------------------------
    diff = mean(AI) - mean(PI)                                    t =  -1.4413
H0: diff = 0                                     Degrees of freedom =      592

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.0750         Pr(|T| > |t|) = 0.1500          Pr(T > t) = 0.9250

. 
. 
. log close
      name:  <unnamed>
       log:  C:\Users\58592ggr\Arbeit\Projects\Active\VotingIncentives\Data&Ana
> lysis\StataFiles\incetives_in_voting_surveys.smcl
  log type:  smcl
 closed on:   4 Jul 2023, 09:26:16
-------------------------------------------------------------------------------
